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Abstract 

Using the non-perturbative renormalization technique, we calculate the renormalization factors 

in 

for quark bilinear operators made of overlap fermions on the lattice. The background gauge field 
[ — ■ is generated by the JLQCD and TWQCD collaborations including dynamical effects of two or 2+1 

o; , 

ON , flavors of light quarks on a 16 x32 or 16 x48 lattice at lattice spacing around 0.1 fm. By reducing 



o 



c3 



the quark mass close to the chiral limit, where the finite volume system enters the so-called e-regime, 



^ ■ the unwanted effect of spontaneous chiral symmetry breaking on the renormalization factors is 

H \ 

suppressed. On the lattices in the conventional p-regime, this effect is precisely subtracted by 
separately calculating the contributions from the chiral condensate. 
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I. INTRODUCTION 



For lattice calculations of operator matrix elements including those of electroweak effec- 
tive Hamiltonian, the operator matching is a necessary step to absorb the difference of the 
renormalization scheme from the conventional continuum one, such as the MS scheme. This 
is necessary for most composite operators except for those protected by some symmetry, e.g. 
the conserved vector current, since the operators are defined with a given lattice action and 
in general divergent in the continuum limit. This operator matching can be done pertur- 
batively and has been done often at the one-loop level, which induces a potential source of 
large systematic error. Given that the strong coupling constant a s is in the range 0.2-0.3, 
a typical size of the two-loop correction is 4-10%. Non-perturbative technique to calculate 
this operator matching is therefore highly desirable to achieve precise calculation of physical 
quantities. 

The Non-Perturbative Renormalization (NPR) method uses the RI/MOM scheme [l| in 
an intermediate step. This scheme is defined for the amputated Green's function in the 
Landau-gauge with an off-shell momentum, which is space-like. Since the matching between 
the RI/MOM and the MS schemes are known to two-loop order in many important operators, 
the method provides a better matching scheme as a whole, though not the entire steps are 
non-perturbative. Moreover, since the perturbative series is in general more convergent in 
the continuum schemes, the remaining uncertainty can be made small to a few percent level. 

Since the method still requires perturbative expansion, the renormalization condition has 
to be applied in the region where non-perturbative effects are sufficiently small. On the other 
hand, one has to avoid large discretization effects that may arise when the renormalization 
scale is too high. Therefore, the renormalization scale fi must satisfy the condition Aqcd ^ 
ji <C Ti /a, where Aq C d stands for the QCD scale and a is the lattice spacing. This region is 
often called the NPR window. 

The non-perturbative effect may be enhanced when the spontaneous chiral symmetry 
breaking (SCSB) occurs and (almost) massless pions arise jj. The reason is that the pion- 
pole contribution in the pseudoscalar channel diverges towards the massless limit and makes 
it difficult to find the NPR window. With the Wilson-type fermions, the problem is severer 
because the error starts at 0(a) and thus the possible window is narrower in the high 
momentum regime. Even with the on-shell (9(a)-improved Wilson fermion, the problem 
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remains since the off-shell amplitude is considered in NPR. With the chirally symmetric 
lattice actions, such as the domain-wall and overlap fermion formulations, the problem 
becomes more tractable because the 0(a) error is absent even in off-shell amplitudes. 

So far, there have been a number of works that calculate the non-perturbative renormal- 
ization factors with the RI/MOM scheme for the domain- wall [2, ^] and for the quenched 
overlap fermions 4], [|], III . 

In this work, we study the non-perturbative renormalization factors with the RI/MOM 
scheme for the quark bilinear operators in unquenched QCD with overlap fermions. Our 
motivation is two-fold. The first is to "provide the renormalization factors corresponding 

n n n 

to the two-flavor [7[ and 2+1-flavor [8|, |9| gauge configurations generated in the large-scale 
dynamical overlap project by the JLQCD and TWQCD collaborations, including the guar 
mass renormalization factor Z m that has been already used in a series of publications 
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12, Q, 14 , H, Q . The second is to study the pion-pole contribution appearing in the 

NPR calculation in detail and demonstrate a method to control the pion-pole effect in a 
reliable manner. 

Since the low-lying eigenmodes of the Dirac operator are expected to dominate the pion- 
pole contribution, it is possible to trace its effect as a function of quark mass by explicitly 
constructing the relevant piece from the low- mode eigenvalues. To be explicit, the pion- 
pole contribution of the form (qq)/p 2 in the operator product expansion contains the chiral 
condensate (qq), which is finite in the vacuum of spontaneously broken chiral symmetry. 
On the lattice of finite volume V, it quickly vanishes as quark mass becomes smaller than 
~ l/T,V, where £ is the chiral condensate in the infinite volume limit. We identify this 
term by explicitly comparing the lattice data of the (inverse) quark propagator with the 
condensate (qq) constructed from the eigenvalues. Thus, this unnecessary term for NPR 
can be identified and subtracted. It means that the pion-pole contribution is no longer a 
problem for the NPR calculation. Clearly, this is possible only when the chiral symmetry is 
preserved on the lattice. Otherwise, the chiral condensate has a bad cubic divergence even 
in the massless limit, hence the identification of its physical contribution is not feasible. 

This paper is organized as follows. We describe the profile of the gauge configurations 
used in this work in Section [Til In Section UTTl we discuss the NPR method and its relation to 
spontaneous chiral symmetry breaking and present our analysis. Results of the calculation 
are given in Section IIII Dl where we summarize all results of the renormalization factor 



3 



available from simple bilinear operators, namely those for the quark mass, the scalar current, 
the tensor operator and the quark field. (The vector and axial vector currents are treated 
independently.) Our conclusion is given in Section HVl 



II. GAUGE CONFIGURATIONS 



In order to make this paper self-contained, we briefly describe the generation of the gauge 



configurations used in this work. We refer [7. 8, 



We use the overlap fermion formulation 
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9] for more complete description. 

1811 on the lattice for both sea and valence 



quarks. The massless overlap-Dirac operator is defined as 



D ov (0) = m (1 + 75 • sgn [H w (-m )}) 



(HI) 



where Hw(-mo) = 75-Diy(— rrio) is the hermitian Wilson-Dirac operator with a large nega- 
tive mass —mo. The massive operator with a bare mass m is constructed from this as 



D n Jm) 



in 



2mo 



D ov (0)+m. 



(11.2) 



We use the Hybrid Monte Carlo (HMC) algorithm 19| to incorporate the fermionic deter- 
minant det[Z) ov (m sea )] (for each flavor) in the path integral. 

Since the overlap-Dirac operator contains the sign function, the corresponding determi- 
nant changes discontinuously on the border of the global topological charge of the gauge 
field configuration, which makes the simulation time-consuming. In order to avoid touching 
the border, where the sign of the lowest eigenvalue of Hw changes, we introduce two extra 
flavors of heavy Wilson fermions such that they produce a factor 

#w(~ m o) 



det 



(11.3) 



-m ) + fi z 

in the Boltzmann weight. Associated (twisted- mass) bosons are also introduced with a 
twisted mass fi. They play a role to minimize the change of the effective gauge coupling 
induced by those extra fermions. Throughout this paper, we choose m = 1.60 and \x = 0.20 
in the lattice unit. As a result, the topological charge Q of the generated gauge configurations 
is fixed to its initial value {^(J. In this work, we choose Q = 0. Although the correct sampling 
of the 6*-vacuum of QCD is spoiled due to the fixed topology, the difference is suppressed 
for large four-volume V, and it is indeed possible to reconstruct the ^-vacuum physics from 



ensemble NF2e NF2p NF3p-a NF3p-b 

Nf 2 2 2+1 

/9 2.35 2.30 2.30 

a' 1 [GeV] 1.776(38) 1.667(17) 1.833(12) 

lattice size 16 3 x 32 16 3 x 32 16 3 x 48 



(m ud ) 0.002 0.015, 0.025, 0.035, 0.015, 0.025, 0.035, 0.015, 0.025, 0.035, 

0.050, 0.070, 0.100 0.050, 0.080 0.050, 0.100 

m s oo oo 0.080 0.100 

0.002, 

m q 0.015, 0.025, 0.035, 0.015, 0.025, 0.035, 0.015, 0.025, 0.035, 0.015, 0.025, 0.035, 

0.050, 0.070, 0.100 0.050, 0.070, 0.100 0.050, 0.080 0.050, 0.100 

#trajectories 2,000 10,000 2,500 2,500 

#step traj. (NPR) 10 100 10 

#step traj. (WTI) see text 20 5 

#low-modes 50 x 2 50 x 2 80 x 2 

# of pi att (Cpiatt) 2 ) 1,375 (30) 1,375 (30) 1,875 (53) 

Relevant papers [10, 11, 12] [13, 15, 22] [16] 

TABLE I: Parameter set for each gauge ensemble NF2e, NF2p, NF3p-a and NF3p-b. Number of 
trajectories are common for all sea quark masses in each ensemble. 



those evaluated by the path integral in a fixed topology [21] . In any case, such finite volume 
effects are irrelevant for the calculation of the renormalization constants considered in this 
work, as it mainly uses the high momentum regime. 

In Table [U we list the parameter set for each gauge ensemble on which we calculate the 
renormalization factors in this work. We performed two-flavor (Nf = 2) and 2+1-flavor 
(Nf = 2 + 1) runs. One of the simulations "NF2e" is in the so-called e-regime of the chiral 
perturbation theory, which corresponds to a very small sea quark mass so that the pion's 
Compton wave length is longer than the lattice extent. The sea quark mass m sea = 0.002 
roughly corresponds to 3 MeV in the physical unit. Other runs at Nf = 2, "NF2p" , are in the 
conventional p-regime, where we take six values of m sea . The 2+1-flavor runs are performed 
at two different values of the strange quark mass, m s = 0.080 ("NF3p-a") and 0.100 ("NF3p- 
b"), so that we can interpolate (or extrapolate) the data to the physical strange quark mass 



5 



afterwards. For each m s , we take five values of sea quark mass corresponding to the up and 
down quarks m u d- 

We employ the Iwasaki gauge action for the gauge part of the lattice formulation. The 
parameter (3 in the action controls the lattice spacing a; we determine the value of the lattice 
spacing from the Sommer scale r by taking r = 0.49 fm as an input after extrapolating 
the lattice data to the chiral limit m sea = or m u d = at a fixed (3. The spatial lattice 
size is 16 3 and the temporal size is 32 and 48 for the two-flavor and the 2+1-flavor runs, 
respectively. 

The valence quark propagator on each ensemble is computed using the multi-shift solver 
at various valence quark masses. For each ensemble in the p-regime, we take the same set 
of masses for the valence quark as that for the sea quarks as listed in Table [H For NF2e, we 
take seven values of valence quark mass: m q = 0.002, 0.025, 0.015, 0.035, 0.050, 0.070, and 
0.100. 

On each gauge configuration fixed to the Landau gauge, we compute the quark propagator 
S ov (x\x src ) = [-D ou (m 9 ) _1 ] X)a ; SI . c , where the location of the source x src is typically fixed at 
the origin. To calculate the renormalization factors, we work in the (four-dimensional) 
momentum space, 

Sovfriatt) = Yl e-^"-^ov(x|x src ). (II.4) 

X 

To avoid possible large discretization error, we restrict the lattice momentum p\ att such that 
its each element py att = 2irn IA /L fJ , does not exceed unity. The numbers of lattice momenta 
satisfying this condition are listed in Table [B Some of them are degenerate in their magnitude 
(Piatt) 2 ; the number of available data point in (pi att ) 2 is also listed in parentheses. On a 
16 3 x 32 lattice, for instance, we have 1,375 different four-momentum from the condition 
—2 < rii < 2 (i = 1,2,3) and —5 < n 4 < 5, and there are 30 different values of (pi a tt) 2 - 
When analyzing the lattice data, we first average over different four-momenta giving an 
identical (piatt) 2 - 
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III. RI/MOM RENORMALIZATION ON THE LATTICE 

A. Renormalization condition and axial- Ward- Takahashi Identity 

We consider flavor non-singlet bilinear operators of the form qTq' with T = 7^, 7^75, /, 
75 and 7^7;,, that we call V, A, S, P and T, respectively. In the following, we may omit the 
prime in q' that indicates that the quark flavor is different from q, but the flavor non-singlet 
operator is always assumed. 

With the exact chiral symmetry of the overlap fermion, these operators are multiplica- 
tively renormalized as 

(qTq) R ( f ,) = Z r (fxa)(qTq) , (III.5) 

where superscripts R and represent the renormalized and bare operators, respectively. 
For divergent operators S, P and T, the renormalized operator may have a dependence on 
the renormalization scale fi. The multiplicative renormalization factor Zr{f^a) then depends 
on the scale fi, too. For the vector and axial-vector currents, the renormalization scale 
dependence is absent because of the current conservation. In the following notation, we 
may drop the dependence on /1, assuming it implicitly. The quark field q is renormalized as 
q R = Z^i^q . 

In the RI/MOM scheme jlj], the renormalization condition is imposed on the amputated 
Green's function Ar(p) = ^Tr[(S '(p))~ Gr(p) (^(p)) -1 T] to satisfy 

A*(p) = Z^)Z r A° r (p) = 1 (III.6) 

at a space-like off-shell momentum p 2 = fi 2 in the chiral limit. Here, the Green's function 
Gr(p) = {lip^Q^l'Wip)) is amputated by the vacuum expectation value of the quark prop- 
agator (S(p)) and projected with an appropriate gamma matrix T. (The 'Tr' denotes the 
trace over the color and spinor indices.) The RI/MOM scheme is defined for the momentum 
configuration that the in-coming and out-going quark momenta are the same p. Since the 
definition involves the external quark field, which is not gauge invariant, the renormalization 
condition depends on the gauge. In the RI/MOM scheme, the Landau gauge is chosen. 

In the RI/MOM scheme, the wave function renormalization Z q is fixed by imposing the 
condition 

12i 



d(s R ( P )] 



d(s(p)y 



(III.7) 
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NF2p 


NF3p-a 


NF3p-b 


^'sea 


_gWTI 




m ud Z^^ 


0.015 


1.37867(61) 


0.015 1.38934(49) 


0.015 1.38968(47) 


0.025 


1.37703(45) 


0.025 1.38709(40) 


0.025 1.38700(36) 


0.035 


1.37412(40) 


0.035 1.38431(32) 


0.035 1.38408(32) 


0.050 


1.37032(33) 


0.050 1.38031(27) 


0.050 1.38019(31) 


0.070 


1.36441(31) 


0.080 1.37196(21) 


0.100 1.36658(26) 


0.100 


1.35436(29) 






0.00 


1.38222(82) 


"chiral limit": 


1.39360(48) 


(xVdof 


0.43 


0.16) 





TABLE II: Summary of the results of Zj TI as a function of m sea or m uc [ for NF2p, NF3p-a and 
NF3p-b. The last two rows show the results of the extrapolation to the chiral limit as described 
in the text. 

at p 2 = fi 2 in the chiral limit. Numerically, though, this is not straightforward since it 
involves a numerical derivative in terms of p^. Instead, we obtain Z q using (IIII.6P for the 
axial- vector vertex function Aa{p) with an input of Za obtained through the axial- Ward- 
Takahashi identity 

Zj T1 A^(A 4 {x)O{0)) = 2m q (P(x)O(0)), (III.8) 

where A4 and P are the axial-vector current in the time direction and pseudo-scalar density, 
respectively. A 4 denotes the symmetrized difference. This relation must be satisfied as far 
as the position x of the operator is not too close to the origin, where some interpolating 
field O is set. Once Za is fixed from this relation, the wave function renormalization is 
determined as Z^(fi) = Z^ ti Aa(p) at p 2 = fi 2 . 

In practice, we use a pseudo-scalar density with a smeared operator for O and sum over 
spatial lattice sites. Then, we fit a ratio 2m q ^2^{P(x,t)O(0)) / A t ^2^{A^(x,t)O(0)) with 
time slices t > t , which is large enough to obtain a constant Z^ T1 . For NF2e, setting 
m q = m seSu = 0.002 and to = 4, we obtain 

Z^ Tl = 1.3511(12). (III.9) 

In other ensembles, Z)^ TI is obtained for each sea quark mass with the valence quark mass 
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sT 137 - 



1.35 



1.34 





FIG. 1: Chiral extrapolation of Z^ Tl for NF2p (left panel) and NF3p-a and NF3p-b (right panel). 



equal to the sea (up and down) quark mass. Results with to = 7 for all ensembles are 
summarized in Table [ITJ where the second row from the last lists the values extrapolated to 
the chiral limit. In the chiral extrapolation, we assume linear plus quadratic dependence on 
m q . Since the local axial- vector current we use on the lattice is not a conserved current at 
finite lattice spacings, the Ward-Takahashi identity ( 1111.81) may be slightly violated. To be 
explicit, a discretization effect of the form a 2 m q d^P is possible as an additive correction to 
A^, which leads to the linear dependence on m q . Including possible quadratic quark mass 
dependence, we use 



zT\m q ) = zr\o) + c imq + CW r 



(in.io) 



by setting the valence quark mass as m q = m sea for NF2p and as m q = m u d for the combined 
data of NF3p-a and NF3p-b. For the case of Nf = 2 + 1, we assume independence of on 
m s , which appears only quark. This assumption is indeed supported by the lattice 

data at two different m s . 

The vertex function Arte) is calculated on the lattice at many different momentum values 
Piatt, whose number is listed in Table [H With the overlap fermion, we compute the vertex 
functions as 



Arte 



latty 
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-Tr 



(-S'ovtelatt)) 1 (<5ov(Platt)r7 5 S'i v (piatt)75) (75^(^)75) ^ (III. 11) 



where the quark propagator is effectively given as 

2rriQ 



S ov (p) 



2m — m q 



-ip-x Brc 



S ov (p) - 



2m n 



(in.i2) 



A. (m =0.0(12) 



9 A 1 t| 



1.15 ■ 



1*4 



T 



I 



• A, 



t 



I I Si in 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



i 



i * 



.15- ^ 



I I 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Pint Ph., 

FIG. 2: Vertex functions A|f t1: (j5 latt ) (circles) and Ay^ipi^t) (squares) at m q = 0.002 for NF2e 
(upper left) and those in the limit of massless valence quark {m q = 0) with the sea quark mass 
fixed at rn se& {m u d) = 0.015 for NF2p (upper right), NF3p-a (lower left) and NF3p-b (lower right). 

This modification of the quark propagator from S ov (p), the inverse of the overlap operator 



D ov (m q ), is made in order to incorporate the quark field rotation q 



'1 - 



£>ov(0) 

2mo 



q,q^q, 



which is necessary to remove the 0(a) effects from off-shell quantities. In ( 1III.11I) . we note 
that 75>S'ov(j'iatt)75 cannot be simply replaced by S'ov(piatt) since the l.h.s of ( 111.41) still depends 
on the source point x src . 



B. Vector and axial-vector vertex functions 



Results for the vector and axial- vector vertex functions are shown in Figure [2] as a function 
of pf att = (ji,a) 2 . In the figure, panels from NF2p, NF3p-a and NF3p-b show the data from 



10 



the lightest m sea or m u d- The chiral symmetry implies that these two functions become 
identical in the massless limit unless the symmetry is spontaneously broken. With exact 
chiral symmetry of the overlap fermion, this should be the case even at finite lattice spacings. 
The result in the e-regime (NF2e, upper-left in the figure) clearly shows this behavior, which 
is consistent with the absence of spontaneous symmetry breaking on a finite volume lattice. 

Other three panels, that are obtained in the p-regime, show the splitting between the 
vector and the axial-vector channels. The numerical data in these plots are naively extrap- 
olated to the chiral limit of the valence quarks by assuming a linear dependence on m q , but 
the qualitative picture remains unchanged for each valence quark mass. 

This inconsistency among the vector and axial-vector currents may be explained as an 
effect of the spontaneously broken chiral symmetry. Even on a finite volume lattice, the 
spontaneous symmetry breaking induces non-zero value of the chiral condensate — (qq) = £ 
as far as the quark mass is much larger than a typical scale 1/EV. An Operator Product 
Expansion (OPE) analysis 3]] suggests that there are contributions of the form Aq CD /p 2 
and mAqcn/p 2 to the difference Aa(p) — Ay(p). These contributions are induced when 
the momentum assignment for the three-point function gives vanishing momentum trans- 
fer at the vertex. Namely, when the incoming and outgoing momenta are identical as in 
the RI/MOM-scheme momentum set-up, which is called the "exceptional momenta", the 
higher dimensional terms in OPE like (qq) 2 /p 6 (with some gamma matrices inserted in the 
numerator) may lead to a much larger contribution of the form ((qq) 2 ) / (Aq CB p 2 ) , which 
remains in the chiral limit in contrast to the lower order contributions m 2 /p 2 or m(qq)/p 4 
[sj]. This problem can be avoided by choosing other momentum configurations, such as the 
RI/SMOM scheme considered in 23]. 

We do not go into details of this problem. But, since the effect becomes statistically 
significant only below pj 2 att ~ 1.0-1.5 for the vector and axial- vector channels, we simply use 
the region that is not largely affected by this effect in the following analysis. 

The quark field renormalization factor Z q (\x) can be obtained from Aa{p) by multiplying 
the axial-current renormalization constant Z^ T1 as determined from the Ward-Takahashi 
identity. The results are shown in Figure [3] by filled circles as a function of p 2 atv The different 
panels represent the data from the ensembles NF2e, NF2p, NF3p-a, and NF3p-b, respec- 
tively. By multiplying the matching factor l/w^ MOM (q) at the four-loop level as defined 
in Appendix [A], we may define the Renormalization Group Invariant (RGI) quantity, which 
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1.6 - 




FIG. 3: Quark field renormalization factor Z q as a function of pjl tt - For NF2p, NF3p-a and 
NF3p-b, data at m sea (m U( i) = 0.015 are plotted as an example. Results in the RI/MOM scheme 
are shown by circles, while those in RGI are plotted by squares. 

is also scheme independent. Our numerical results plotted by squares in Figure [3] clearly 
show the expected scale independence. Since we expect discretization effects proportional 
to a 2 £>Ltt' we n t the lattice data above (pi a tt) 2 = 1-0 by a linear function and obtain the 
result for Z^ Gl from an intercept at p 2 att = 0. The lattice data below a 2 Pi aU — 1.0 are largely 
affected by the effect of spontaneously broken chiral symmetry and deviate from the linear 
behavior as expected. 

The results for Z^- G1 are listed in Table IIHI Also listed are the results converted to the 
MS scheme at \x = 2 GeV using the four-loop level matching constant u^ s (/i) defined in 
Appendix [A] 
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0.00 0.01 0.01 0.01 0.02 0.03 0.03 0.04 
2m . +m 



FIG. 4: Chiral extrapolation of Zf s {2 GeV) on the NF2p (left) and NF3p-a (right, filled symbols) 
and NF3p-b (right, open symbols) lattices. 

So far, the results are given at each sea quark mass after taking the chiral limit of valence 
quarks. The chiral limit of sea quarks can be taken by assuming that the sea quark mass 
dependence has the form Z(l + 2c (2) m^ ca ) (for NF2p) or Z(l + c (3) (2m^ + m 2 s )) (for NF3p-a 
and NF3p-b). The coefficients and are numerical constants depending on the number 
of flavors. The linear term in m sea (or in my) should not remain for the quantities irrelevant 
to the chiral symmetry breaking. Figure H] shows the chiral extrapolation of Z q ^ ls (2 GeV) 
for both NF2p and NF3p-a/NF3p-b. We do not observe any significant sea quark mass 
dependence. The chiral extrapolation should therefore be very stable. The results are listed 
in Table HH 



C. Scalar and Pseudo-scalar vertex functions 

In Figure the momentum dependence of the scalar vertex function A l g U (pi a tt) (filled 
symbols) and the pseudo-scalar vertex function Ap tt (pi at t) (open symbols) is shown for each 
ensemble. For the data in the e-regime (circles in the upper left panel), we observe an 
excellent agreement between As(p) (filled symbols) and Ap(p) (open symbols), which is 
expected from the exact chiral symmetry of the overlap fermion. On the other hand, once 
the valence quark mass m q is out of the e-regime (the data at m q = 0.015, 0.025 and 0.035 
are plotted by squares, diamonds and triangles, respectively), we find large disagreement 
between As(p) and Ap(p). 

This observation again indicates the effect of the spontaneous chiral symmetry breaking. 
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m sea (m ud ) 


Z HGI 


Z g MS (2 GeV) 


NF2e, 0.002 


1.4170(47) 


1.4799(50) 


NF2p, 0.015 


1.4540(54) 


1.5186(56) 


0.025 


1.4503(51) 


1.5147(53) 


0.035 


1.4479(53) 


1.5122(56) 


0.050 


1.4486(49) 


1.5129(51) 


0.070 


1.4442(51) 


1.5083(54) 


0.100 


1.4500(59) 


1.5143(62) 


0.000 


1.4526(30) 


1.5170(31) 


(x7dof 


0.33 


0.33 ) 


NF3p-a, 0.015 


1.4575(30) 


1.5279(31) 


0.025 


1.4641(38) 


1.5348(40) 


0.035 


1.4660(51) 


1.5368(54) 


0.050 


1.4528(27) 


1.5230(29) 


0.080 


1.4590(37) 


1.5294(39) 


NF3p-b, 0.015 


1.4565(40) 


1.5269(42) 


0.025 


1.4585(32) 


1.5290(34) 


0.035 


1.4555(28) 


1.5258(29) 


0.050 


1.4467(30) 


1.5166(31) 


0.100 


1.4578(55) 


1.5283(57) 


"chiral limit" : 


1.4592(29) 


1.5296(31) 


(x 2 /dof 


2.20 


2.20 ) 



TABLE III: Numerical results for the quark wave function renormalization factor Z q . The values 
in the RGI definition Zf Gl and those defined in the MS scheme at fi = 2 GeV, Zf s (2 GeV) are 
listed for each sea quark mass. 

From the OPE analysis one expects that this effect is more significant than in A v (p) and 
Aa(p), because the violation is enhanced by an inverse quark mass as discussed below. 
One has to subtract this effect in order to extract the renormalization constants because 
its matching is based on the continuum perturbation theory that does not contain non- 
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FIG. 5: Vertex functions A l g U (pi att ) (filled symbols) and Ap tt (pi att ) (open symbols). The data are 
shown for NF2e (upper left), NF2p (upper right), NF3p-a (lower left) and NF3p-b (lower right) 
as functions of pf att - In each panel, data at m q = 0.015, 0.025 and 0.035 are shown by squares, 
diamonds and triangles, respectively. For NF2e, data at m q = 0.002 are shown as well by circles. 
For NF2p, NF3p-a and NF3p-b, data at the lightest sea quark m seSb (m uc i) = 0.015 are plotted as 
an example. 



perturbative effects. 

We consider the quark mass dependence of As(p) and Ap(p) using OPE along the line 
of the analysis in . Using the vector and axial- vector Ward-Takahashi identities, one may 
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obtain relations between the vertex functions and the inverse quark propagators as 24] 

On the lattice we use the improved overlap quark propagator S nv (p) in place of From 
OPE the inverse quark propagator Tr (S(p)) _1 may be written as 25] 



-3-Tr (Sip))- 1 = C • M + Z g Z m m, + • ■ ■ (111.15) 
12 £r 



in the large p 2 regime. The effect of the chiral symmetry breaking is picked up through the 
chiral condensate (qq), and C is a perturbatively calculable constant. At the one-loop level, 
C = 47ra s /3. As in the case of the vector and axial- vector vertex functions, the effects from 
higher dimensional operators, such as ((qq) 2 ), may also exist. They are usually suppressed 
by additional powers of 1/p 2 , but due to the lack of the momentum injection the suppression 
may not work in the case of the inverse quark propagator and of the vertex functions at 
zero momentum transfer. We therefore leave C as an unknown constant instead of using the 
perturbatively known value. 

Then, using the relations ( 1III.13I) and ( 1III.14I) . we may evaluate the effect of the chiral 
condensate on the vertex functions as 



A s (p) = + Z q Z m + ■■■ , (111.16) 

Ap(p) = - 2 ^ + Z q Z m + --- . (111.17) 

p z m q 

From these expressions, one sees an enhancement in the low p 2 region due to the chiral 
condensate only for the pseudo-scalar channel, while the scalar channel should not be affected 
too much because of a derivative with respect to m q rather than a factor l/m q . 

The quark mass dependence of the condensate (qq) is not a trivial issue, since it has 
effects from both ultraviolet and infrared origins. Since the operator qq contains quadratic 
divergence of the form m q /a 2 apart from the chiral limit, the chiral condensate (qq) directly 
calculated on the lattice contains unphysical large m q dependence. It has to be subtracted 
before the analysis using flIII.16j) and (IIII.17I) . because the formulae are obtained as an 
expansion around the chiral limit. 
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In the infrared regime, the chiral condensate has a non-trivial quark mass dependence 
especially in a finite volume. First, because of the pion-loop effects, the chiral condensate 
develops the chiral logarithm of the form m q In m q with known coefficients 26( . On a finite 
volume lattice, the quark mass dependence becomes more complicated. Namely, once the 
quark mass enters the e-regime, the mass dependence is no long er g overned by the simple 
chiral logarithm, but given by the formula recently developed in 271 ]. 

In our analysis, instead of using the formula in 27] we calculate the condensate using 
its eigenvalue decomposition by making use of the low eigenmodes obtained on the same 
ensembles. For each lattice configuration, we define 



(qq) 



N 



2m n 



LIU j^ml + A* A, 



(111.18) 



where A« is an eigenvalue of the massless overlap-Dirac operator, which satisfies the eigen 
equation 



Aw(0) 1 



£>ov(0) 
2mn 



-i 



Ui(x) = XiUi(x) 



(111.19) 



with Ui(x) an eigenvector. In (IIII.18I) we use the fact that the eigenvalues appear as complex 
conjugate pairs. The normalization in ( 1111.181) contains the lattice volume L^Lf. 

We truncate the sum in ( 1III.18I) at iV-th eigenvalue, which may be considered as a "renor- 
malization scheme" to define the divergent operator qq. Here, N plays a role of the ultraviolet 
cut-off. After taking an ensemble average, we denote the chiral condensate thus defined as 
(qq)( N \ In the course of our project, we calculate and store the low-lying eigenvalues and 
eigenvectors of the overlap-Dirac operator. In addition to the calculation of the truncated 
chiral condensate ( 1111.181) . these eigenmodes can be used to precondition the solvers, to aver- 



age over source 
observables 
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joints, or to construct disconnected diagrams in the calculations of physical 
. The numbers of the stored low-modes for each configuration are 
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listed in Table [fl 

From a dimensional analysis, the quark mass dependence of (qq) ( N > may be parametrized 

as 

(qq)W = (qq)^ ht) + cf ^ + cf^ . (111.20) 

Because of the exact chiral symmetry of the overlap fermion, there is no leading power 
divergence of the order 1/a 3 , and the term behaves as m 2 q /a is also absent. Although the 
cubic term c^m q in (lIII.20j) may accompany a logarithm lnm ? , we omit it for simplicity 
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FIG. 6: Comparison between (qq) obtained from all-to-all propagator (circles) and (qq)^ N 50 ^ 
(squares) as a function of m sea = fn q . The data obtained on NF2p are shown. 

as the m q term itself is a minor correction. The subtracted condensate (gg)( subt ) is then free 
from power divergences, but could still contain non-divergent m q dependence, such as the 
chiral logarithm. 

In FigureEl we compare a "full" calculation of (qq) (circles) corresponding to iV = Y2L\L t 
and (qq)^ N ~ 50 ^ ■ The "full" calculation contains the contributions of all eigenmodes which are 



evaluated by a stochastic method. (Our set-up is explained in 28].) The data on the NF2p 
lattice with sea and valence quark masses set equal are shown as an example. The results 
clearly show that the divergent term m q /a 2 in (IHI.20I) dominates the "full" condensate, and 
it seems difficult to extract (qq)( subt ^ from this data alone. The truncated condensate, on 
the other hand, does not have that strong m q dependence, but still both the m q /a 2 and m q 
terms are visible. 

We now try to extract the non-divergent term (gg)( subt ) using (lIII.20j) . In Figures [7H9] 
(left panel) we plot the truncated condensate (qq) as a function of the valence quark mass 
with three or four different values of N. The data are shown for individual lattice ensembles 
(NF2e, NF2p, NF3p-a and NF3p-b); except for NF2e the results at the lightest sea quark are 
shown as an example. The truncated condensate can be constructed at arbitrary values of 
the valence quark mass m q without extra computational costs. In order to see the ultraviolet 
behavior, we plot in the mass region up to m q = 0.20, which is twice larger than the largest 
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FIG. 7: Subtraction of the power divergence in the chiral condensate. The left panel shows 
(qq)^ calculated on the NF2e lattice as a function of the valence quark mass m q (solid curves). 
The number of low-modes included are N = 50, 40 and 30 from top to bottom. The solid curves 
are data constructed from calculated eigenvalues, and the points with error bars are representative 
points used in our fit. The fit curves according to (|III.20j) are shown by dashed curves. The right 
panel represents the subtracted condensate (qqp anbt \ 
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q 4 

FIG. 8: Same as Figure but for NF2p at m sea = 0.015. 

simulated sea quark mass. 

When we fit the lattice data to (jHL20p . we take five or six representative points of m q in 
the range 0.10 < m q < 0.18 for NF2e and NF2p or 0.10 < m q < 0.20 for NF3p-a and NF3p- 
b. The upper limit is chosen such that |Ajv| > m q for the given N. Otherwise we do not 
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FIG. 9: Same as Figure \7\ but for NF3p-a (top) and NF3p-b (bottom) at m u d = 0.015. The 
number of low-modes included are N = 80, 70, 60 and 50 from top to bottom. 



expect the ultraviolet behavior ( 1III.20I) . In this rather heavy mass region, we do not expect 
additional mass dependence from the infrared origin, and we simply set — (qq)^ snbt ' = S 
with E a constant. 

The fit results are shown in Figures [7H9] (left panel) by dashed curves. In each right panel 



of these figures, the subtracted condensates (qq) 



(subt) 



(qq) 



( N ) J N ) /„2 



c\"' /a 2 m q — c^'nvL for the 



same choices of N as the left panel are shown. We observe that the subtracted condensate 
depends on N very mildly to a few % order. It implies that our subsequent analysis using 
(qq) ^ with the maximal iV may contain a small systematic error due to the truncation of 
N. We discuss this point and estimate the error in Section IIII D[ 

We use (qq)( suht ^ thus obtained at each sea quark mass as a function of the valence 
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FIG. 10: Vertex functions A^ tt (pi a tt) (filled symbols) and Ap (piatt) (open symbols). The data are 
shown for the NF2e (upper left), NF2p (upper right), NF3p-a (lower left) and NF3p-b (lower right) 
lattices as functions of m q . In each panel, data at selected values of p? tt and their fit curves are 
presented. For NF2p, NF3p-a and NF3p-b, data with a fixed sea quark mass m sea (m n d) = 0.015 
are plotted as an example. 

quark mass in the analysis of the scalar and pseudo-scalar vertex functions, (IIII16I) and 
(IIII.17p . respectively. The valence quark mass dependence of A^Qo^) and Ap tt (pi at t) at 
four representative values of pf att are plotted in Figure [TOl We find that both the scalar (filled 
symbols) and pseudo-scalar (open symbols) vertices are nicely described by the fit curves 
according to (1IH.16|) and (1III.17j) supplemented by the measured {qq)^ snht \ In particular, as 
seen near the chiral limit of the NF2e data, the fit curves precisely reproduce the agreement 
of Ag t1: (p latt ) and Ap^Qoiatt) in the e- regime, which is not expected when (qq) is treated as a 
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pj tt C Z q Z rn B s B P X 2 /dof 

1.426 5.92(40) 1.443(16) -3.10(90) -1.35(34) 0.025 

1.889 5.01(32) 1.384(13) -2.08(62) -0.85(25) 0.017 

2.352 4.32(28) 1.342(11) -1.60(48) -0.76(22) 0.016 

2.814 3.78(26) 1.3060(92) -1.23(36) -0.65(15) 0.017 

TABLE IV: Parameters in the simultaneous fit of A 1 ^** and Ap tt for NF2e. Results at representative 
values of lattice momenta are listed. 

mass- independent constant. 

In addition to fjlll. 16|) and 1T[) . quadratic mass-dependence is possible for the vertex 
functions A^ tt (p latt ) and A£ tt (p Ultt ): 

As(piatt) = + Z * Z m + B s m 2 q , (111.21) 

Ap(piatt) = ^^- + Z q Z m + B P m 2 . (111.22) 

From a combined fit of the valence quark mass dependence, we obtain the parameters C, 
Z q Z m , Bp and Bs at each value of pf att . Numerical results are listed in Tables ITVl - IVHI for 
each sea quark masses of the NF2e, NF2p, NF3p-a and NF3p-b lattices. We find that the 
values of C depend on p 2 only mildly, which is consistent with the logarithmic dependence 
through ot s . 

D. Renormalization constants 



From the fits described in the previous subsections, we obtain the numerical results for Z q 
and Z q Z m = Z q /Zs for each available values of pi a tt- From a similar analysis, we also obtain 
At = ZgZ^ 1 , which does not depend on the quark mass significantly. We combine these 
results with Z q (fi) to obtain Zm^ M ° M {fi) = l/Z^' I/ ' MOM (/i) and Z^ MOM (fi) as functions of 
the renormalization scale \x. 

The results are plotted in Figures [11] and [12] for Z m and Z T) respectively. Filled black 
symbols representing the numerical data for the RI/MOM scheme clearly show a scale (or 
Piatt) dependence. This dependence can partly be absorbed by perturbatively calculated 
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2 

m sca, Piatt 


C 


7 7 
^q^m 


B S 


Bp 


X 2 /dof 


0.015 1.426 


8.48(83) 


1.472(22) 


-2.0(1.9) 


-0.80(91) 


0.004 


1.889 


7.16(68) 


1.412(16) 


-1.5(1.3) 


-0.57(66) 


0.005 


2.352 


6.19(59) 


1.366(11) 


-1.01(95) 


-0.35(46) 


0.005 


2.814 


5.31(51) 


1.3271(89) 


-0.65(75) 


-0.18(38) 


0.004 


0.025 1.426 


7.43(39) 


1.497(20) 


-5.3(1.1) 


-2.21(49) 


0.139 


1.889 


6.29(32) 


1.429(15) 


-3.53(76) 


-1.33(35) 


0.084 


2.352 


5.43(26) 


1.380(13) 


-2.66(56) 


-1.07(25) 


0.071 


2.814 


4.61(23) 


1.341(10) 


-2.18(43) 


-0.96(23) 


0.094 


0.035 1.426 


8.57(91) 


1.446(24) 


-0.2(2.0) 


-0.1(1.1) 


0.190 


1.889 


7.37(77) 


1.392(17) 


0.1(1.4) 


0.11(71) 


0.179 


2.352 


6.23(61) 


1.353(13) 


-0.13(92) 


-0.02(48) 


0.133 


2.814 


5.40(56) 


1.317(11) 


0.00(76) 


0.09(43) 


0.129 


0.050 1.426 


9.00(58) 


1.469(21) 


-1.3(1.5) 


-0.86(63) 


0.021 


1.889 


7.61(49) 


1.419(15) 


-1.1(1.1) 


-0.54(48) 


0.014 


2.352 


6.54(42) 


1.377(12) 


-0.95(81) 


-0.59(36) 


0.019 


2.814 


5.60(37) 


1.340(10) 


-0.89(63) 


-0.61(31) 


0.047 


0.070 1.426 


7.33(54) 


1.476(19) 


-5.1(1.4) 


-2.19(63) 


0.303 


1.889 


6.23(45) 


1.412(14) 


-3.56(99) 


-1.50(45) 


0.267 


2.352 


5.36(40) 


1.363(11) 


-2.61(78) 


-1.03(38) 


0.228 


2.814 


4.56(34) 


1.3262(93) 


-2.23(58) 


-0.94(30) 


0.249 


0.100 1.426 


8.24(61) 


1.444(19) 


-1.9(1.5) 


-1.09(62) 


0.008 


1.889 


6.99(53) 


1.394(14) 


-1.4(1.1) 


-0.81(46) 


0.014 


2.352 


6.04(47) 


1.351(12) 


-0.92(88) 


-0.51(42) 


0.003 


2.814 


5.18(42) 


1.316(10) 


-0.72(74) 


-0.32(39) 


0.004 



TABLE V: Same as Table ED but for NF2p. 

matching factor 

w RI/MOM (/i) {<D 

= m or T) to the RGI values as in the case of Z q . The 
perturbative results for w^ MOM (fi) and w^ s (fi) are summarized in Appendix lAl 

The numerical data for Zg GI = Z™ /MOM (p latt )/w™ /MOM (p latt ) are also shown in Fig- 
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m ud Putt 


C 


7 7 
^q^m 


B S 


Bp 


X 2 /dof 


0.015 1.405 


9.00(98) 


1.445(21) 


2.3(2.8) 


1.3(1.5) 


0.345 


1.851 


7.58(82) 


1.393(14) 


1.3(1.9) 


0.9(1.1) 


0.248 


2.279 


6.65(72) 


1.347(10) 


1.2(1.5) 


0.87(86) 


0.276 


2.690 


5.86(62) 


1.3228(83) 


0.8(1.1) 


0.55(58) 


0.194 


0.025 1.405 


7.18(52) 


1.477(14) 


-3.5(1.4) 


-1.59(60) 


0.046 


1.851 


6.14(44) 


1.418(11) 


-2.77(94) 


-1.41(42) 


0.056 


2.279 


5.25(37) 


1.3661(93) 


-2.07(72) 


-1.03(35) 


0.046 


2.690 


4.64(32) 


1.3379(78) 


-1.86(53) 


-0.92(26) 


0.065 


0.035 1.405 


7.21(54) 


1.495(19) 


-3.3(2.0) 


-1.90(81) 


0.018 


1.851 


6.19(46) 


1.427(15) 


-2.1(1.5) 


-1.17(58) 


0.018 


2.279 


5.36(40) 


1.376(12) 


-1.6(1.1) 


-0.91(48) 


0.015 


2.690 


4.75(36) 


1.3448(97) 


-1.20(88) 


-0.58(38) 


0.010 


0.050 1.405 


7.88(56) 


1.444(16) 


-0.2(1.9) 


-0.34(82) 


0.098 


1.851 


6.80(52) 


1.391(12) 


0.0(1.5) 


-0.03(72) 


0.099 


2.279 


5.82(44) 


1.3464(92) 


0.0(1.1) 


0.14(52) 


0.086 


2.690 


5.27(43) 


1.3216(80) 


0.08(98) 


0.10(48) 


0.075 


0.080 1.405 


9.24(67) 


1.449(19) 


-0.7(2.6) 


-9.35(78) 


0.211 


1.851 


8.00(59) 


1.395(13) 


-0.2(1.9) 


-5.81(59) 


0.226 


2.279 


6.81(50) 


1.352(10) 


-0.1(1.4) 


-3.93(47) 


0.189 


2.690 


6.06(44) 


1.3296(86) 


-0.2(1.1) 


-3.01(38) 


0.138 



TABLE VI: Same as Table ED but for NF3p-a. 



ures [TT] and [121 We find that the scale dependence is largely absorbed at least above 
(piatt) 2 — 1, as expected. Below (pi a tt) 2 — 1 the perturbative estimate of w^^ MOM (p) be- 
comes less precise even though three- or four-loop calculations are used. Remaining scale 
dependence above (piatt) 2 — 1 is ascribed to the discretization effect of 0(a 2 ). We therefore 
extrapolate the data for Z@ above (piatt) 2 = 1 to the vanishing (p\ a tt) 2 limit assuming a 
linear dependence on (pi a tt) 2 ; which is shown by solid lines in Figures [TT1 and [T2l 

The renormalization constants in the MS scheme are obtained as Z^ s (/i) = w^ s (fi) -Z^ Gl , 
again using the matching factor to the RGI value WQ S (fi). Results of the RGI value and 
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m ud Putt 


C 


7 7 


B S 


Bp 


X 2 /dof 


0.015 1.405 


8.08(59) 


1.464(13) 


-2.30(94) 


-0.99(55) 


0.020 


1.851 


6.92(49) 


1.4009(98) 


-1.57(63) 


-0.57(35) 


0.029 


2.279 


5.97(42) 


1.3527(78) 


-1.04(50) 


-0.34(29) 


0.010 


2.690 


5.30(37) 


1.3239(68) 


-0.76(41) 


-0.19(22) 


0.007 


0.025 1.405 


7.04(42) 


1.4699(93) 


-2.61(78) 


-1.21(42) 


0.012 


1.851 


6.01(36) 


1.4101(75) 


-1.88(56) 


-0.83(31) 


0.004 


2.279 


5.19(32) 


1.3595(60) 


-1.34(43) 


-0.54(22) 


0.003 


2.690 


4.65(29) 


1.3321(52) 


-1.11(34) 


-0.51(18) 


0.005 


0.035 1.405 


7.14(50) 


1.475(11) 


-3.44(78) 


-1.39(43) 


0.146 


1.851 


6.09(41) 


1.4126(85) 


-2.47(54) 


-1.07(27) 


0.087 


2.279 


5.26(35) 


1.3624(68) 


-1.84(41) 


-0.77(22) 


0.095 


2.690 


4.67(31) 


1.3320(58) 


-1.31(32) 


-0.48(18) 


0.059 


0.050 1.405 


8.93(56) 


1.423(15) 


-0.1(1.3) 


0.03(55) 


0.295 


1.851 


7.59(47) 


1.3755(99) 


-0.09(84) 


0.03(37) 


0.325 


2.279 


6.52(41) 


1.3363(77) 


-0.16(63) 


-0.01(29) 


0.200 


2.690 


5.85(36) 


1.3115(65) 


-0.04(50) 


0.02(22) 


0.240 


0.100 1.405 


8.68(56) 


1.468(16) 


-1.7(1.2) 


-0.61(60) 


0.123 


1.851 


7.34(48) 


1.412(12) 


-1.21(91) 


-0.49(46) 


0.113 


2.279 


6.26(40) 


1.3650(87) 


-1.00(66) 


-0.41(29) 


0.072 


2.690 


5.61(35) 


1.3353(75) 


-0.59(54) 


-0.24(26) 


0.104 



TABLE VII: Same as Table HV] but for NF3p-b. 



those in the MS scheme at [i = 2 GeV are listed in Table IVIIII for O = m and S with the 
four-loop level matching, and in Table HXl for O = T with the three-loop level. 

For the NF2p ensembles the renormalization factors at finite sea quark masses are ex- 
trapolated to the limit of m sea = as a linear function of 2m^ ea . For the 2+1-flavor data, we 
combine NF3p-a and NF3p-b to quote the final result in the chiral limit of all of the three 
flavors, assuming a sea quark mass dependence of the form Z(l + c^\2m? ud + m? s )). The 
extrapolation is shown in Figures [TBI and fT4l for the NF2p (left panel) and NF3p-a/b (right 
panel) ensembles. Although we do not observe any systematic sea quark mass dependence, 
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FIG. 11: Renormalization factors for the quark mass Z m in RI/MOM scheme (circles) and their 
RGI values (squares) as functions of pf att for NF2e (upper left), NF2p (upper right), NF3p-a (lower 
left), and NF3p-b (lower right). Results of the linear extrapolation to the (pi a tt) 2 ~~ * limit 
of are shown as well. For NF2p, NF3p-a and NF3p-b, data with a fixed sea quark mass 

m sea (m u rf) = 0.015 are plotted as an example. 

the data show larger fluctuations than the statistical errors at each sea quark mass for Z m . 
As a result, the x 2 /dof for the combination of NF3p-a and NF3p-b is uncomfortably large 
(~ 2.6), as listed in Tables IVIIIl This may indicate that the statistical error estimated at 
each sea quark mass is underestimated. It is also suggested from the size of the statistical 
error at a fixed sea quark mass, say m sea (or m U( i) = 0.015. Namely, the size of error is 
comparable between NF2p and NF3p-a/b, though the statistics is more than factor of two 
larger for NF2p. We use the jackknife method for the statistical analysis with a bin size of 
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FIG. 12: Same as Figure [TT] but for the tensor current renormalization factor Zt- 
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FIG. 13: Left panel: Z^ s (2 GeV) for NF2p as a function of 2™^. Right panel: Same value for 
NF3p-a and NF3p-b as a function of 2m^ rf + ml . Linear extrapolation to the m sea = limit or the 
m-tid = limit is shown as well. 
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m sea ,(m ud ) 


zf GI 


vRGI 


Zf s (2 GeV) Z™ s (2 GeV) 


NF2e, 0.002 


0.709(11) 


1.411(21) 


1.205(18) 


0.830(12) 


NF2p, 0.015 


0.743(18) 


1.345(33) 


1 263(30) 


791(19) 


0.025 


0.719(12) 


1.390(24) 


1 99^91 ^ 




0.035 


0.764(20) 


1.308(37) 


1.298(35) 


0.769(22) 


0.050 


0.746(16) 


1.339(29) 


1 268(27) 


788(17) 

KJ • 1 K^JK^J \ -L- 1 / 


0.070 


0.726(14) 


1.378(26) 


1.234(23) 


0.810(15) 


0.100 


0.761(14) 


1.313(26) 


1.293(24) 


0.772(15) 


0.000 


0.7309(87) 


1.366(16) 


1.243(15) 


0.8035(97) 


(xVdof 


1.15 


1.15 


1.15 


1.15) 


NF3p-a, 0.015 


0.766(19) 


1.303(34) 


i.zyo^oz ) 


f) 770/ / 9rA 


0.025 


0.734(10) 


1.362(20) 


1.242(17) 


0.805(12) 


0.035 


0.721(17) 


1.386(33) 


1.221(29) 


0.819(19) 


0.050 


0.757(15) 


1.320(27) 


1 281(25) 

J- • — KJ J- \ — '.J 1 


780(16) 

KJ • 1 K^J KJ V 1 V/ f 


0.080 


0.765(17) 


1.304(31) 


1 296(29) 


770(18) 

W • 1 1 KJ \ i- W / 


NF3p-b, 0.015 


0.748(11) 


1.337(19) 






0.025 


0.7354(62) 


1.359(12) 


1.245(11) 


0.8031(70) 


0.035 


0.7311(87) 


1.368(17) 


1.238(15) 


0.8078(98) 


0.050 


0.774(13) 


1.289(23) 


1.310(22) 


0.761(14) 


0.100 


0.748(13) 


1.336(25) 


1.266(22) 


0.789(15) 


"chiral limit" : 


0.7325(88) 


1.364(16) 


1.240(15) 


0.8057(97) 


(xVdof 


1.61 


1.64 


1.61 


1.64) 



TABLE VIII: Renormalization factors of the scalar operator and quark mass in the RGI and in 
the MS schemes at jj, = 2 GeV. The results at each sea quark mass are listed as well as those in 
the chiral limit of sea quarks. 

50 HMC trajectories. Given the limited total length of trajectories (2,500 for NF3p-a/b), 
the statistical error does not change much even if we increase the bin size to 100 trajecto- 
ries. We do not investigate this point further, because the statistical error does not give the 
dominant part of the error in the final results. 
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m sca (m ud ) 


7 RGI 

Zjrp 


Z™ s {2 GeV) 


NF2e, 0.002 


1.7023(62) 


1.4689(53) 


NF2p, 0.015 


1.7461(69) 


1.5066(59) 


0.025 


1.7418(63) 


1.5030(54) 


0.035 


1.7393(70) 


1.5008(61) 


0.050 


1.7470(72) 


1.5075(62) 


0.070 


1.7361(63) 


1.4981(54) 


0.100 


1.7330(64) 


1.4953(55) 


0.000 


1.7441(38) 


1.5050(33) 


(x 2 /dof 


0.24 


0.24) 


NF3p-a, 0.015 


1.7662(44) 


1.5283(38) 


0.025 


1.7663(50) 


1.5284(43) 


0.035 


1.7685(54) 


1.5303(47) 


0.050 


1.7591(38) 


1.5222(33) 


0.080 


1.7674(48) 


1.5294(41) 


NF3p-b, 0.015 


1.7620(43) 


1.5247(37) 


0.025 


1.7630(42) 


1.5255(37) 


0.035 


1.7637(42) 


1.5261(36) 


0.050 


1.7518(43) 


1.5159(37) 


0.100 


1.7640(56) 


1.5264(48) 


"chiral limit": 


1.7639(35) 


1.5262(30) 


(xVdof 


1.18 


1.18) 



TABLE IX: Same as Table Em] but for Z T . 



For the central values of the final results, we quote the result at m sea = 0.002 for NF2e 
and that in the m sea = limit for NF2p or in the m U( i = m s = limit for the combination 
of NF3p-a and NF3p-b. In Table IVIIIl the extrapolated values are listed in separated rows. 
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FIG. 14: Same as Figure [TJ but for Z| IS (2 GeV). 



The results for Z m are 



Z^{2 GeV) = < 



0.824(14) (24) for iV> = 2, /? = 2.35 

0.804(10) (25) for JV) = 2, = 2.30 

0.806(12)(24) (+°°) for JV/ =2 + 1,^ = 2.30, 



(111.23) 



The first error is statistical, which includes the small statistical errors in the extraction of 
Z^ T1 and the lattice scale a -1 . The scale affects the determination of the matching point 
H — 2 GeV. On this error, we also take account of the ambiguity in removing the scale 
dependence of Z^ Gl by comparing the results with different ranges of the linear fit. The 
systematic errors given in the second and third parentheses are described in the following. 

An important source of the systematic error is the truncation of the perturbative expan- 
sion in the matching between the RI/MOM and MS schemes. It is given by a ratio of two 
matching factors to the Renormalization Group Invariant (RGI) value, i.e. Wm 1 M (/i) and 
w m S (^) i n QA.ip . The perturbative expansion of these factors is given in (1A.2j) and known to 
four-loop order. By setting \i — 2 GeV, we may evaluate how it depends on the loop order. 
For Nf = 2, the ratio io^(//)/wm /MOM (>) becomes 1, 0.911, 0.863, and 0.835 when the 
perturbative expansion includes 0(a°), 0(a]), 0(af), and 0(af) terms respectively. From 
this observation, we find that the perturbative expansion converges such that the additional 
correction is about 60% of the correction of the previous order. Same level of the correction 
is observed for the case of Nf = 2 + 1. We therefore assume that this convergence persists 
at the next unknown perturbative coefficient. The second error in (1III.23I) is estimated by 
taking a difference of the current best four-loop analysis and the second best three-loop 
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analysis and multiplying a factor 0.6. 

The effect of SCSB may arise in two different ways. First, in the extraction of Z q we used 
the axial- vector vertex function A^, but if we used the vector vertex function Ay instead 
the result is slightly shifted, which is given in the third parentheses. Note that this does 
not matter for NF2e, because there is no significant difference between Ay and A^ in the e- 
regime. Second, one may expect some uncertainty in the process of subtraction of the power 
divergent piece from (qq)^. In Section llH C\ we demonstrate that the power-divergent term 
can be removed from (qq) to obtain (qqy snbt ' in almost ^-independent way. However, it 
does not guarantee that the results for (qq)^ snbt ^ are unchanged beyond the maximum value 
of A" we studied. In fact, in the p-regime, we find that E = (qq)^ anht ^ (m q = 0) obtained 



with various values of A" slightly increases as a function o 
the chiral perturbation theory on the same ensembles 10 



N. In the calculation based on 



29j, we find ~ 10% larger values 
of E than those from the eigenvalue decomposition for NF2p and NF3p-a/b. For NF2e, the 
result of the calculation based on the chiral random matrix theory [U|] is ~ 10% smaller. 
To estimate the effect of the truncation of eigenvalues, we repeat the same analysis by fixing 
the value of E 10% smaller (larger) than the original one for NF2e (NF2p and NF3p-a/b). 
As a result, we find the magnitude of finite A" effect is similar to the statistical errors for 
all cases. We quote the difference from the central value in the third error in (1III.23j) . For 



NF2p and NF3p-a/b, we combine this error with the effect from the difference between A^ 
and Ay, which is in the same direction. 

For completeness, we also present results for Zt- Matching procedures are illustrated in 
Figure H21 Table [IX] summarizes the RGI and the MS values. The left and right panels 
in Figure [14] show the linear extrapolation as a function of 2m;? ea (NF2p) and 2m^ d + m 2 s 
(NF3p-a/b), respectively. 

IV. CONCLUSION 

We calculated the renormalization factors for the quark bilinear operators constructed 



from the overlap fermion formulation, based on the original idea of NPR proposed in [I] 
The aim of this calculation is to provide the renormalization factors for a series of numerical 
studies being performed by the JLQCD and TWQCD collaborations using dynamical overlap 
fermions. By virtue of the exact chiral symmetry of the overlap fermion, the analysis is 
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largely simplified compared to other non-chiral fermion formulations. 

Through the simulation in the e-regime, we explicitly confirm that the vector and axial- 
vector vertex functions agree with each other when the effect of spontaneous chiral symmetry 
breaking is negligible. This may provide a clean way to calculate the renormalization factors 
through the NPR method, since the calculation does not suffer from the potential problems 
due to pion poles. 

In the p-regime, where the spontaneous symmetry breaking effectively remains even on a 
finite volume lattice, we may precisely control the non-perturbative quark mass dependence 
of the quark propagator and vertex functions using the OPE analysis supplemented by the 
condensate explicitly constructed from the low-lying quark eigenmodes. The exact chiral 
symmetry of the overlap fermion plays an important role also in this analysis. 

Our main results are those of the mass renormalization factor Z m , which is an inverse 
of the scalar density renormalization factor Z$- The result has been already used in the 



calculation of the chiral condensate in two-flavor QCD from the Dirac operator spectrum 11 



and from the topological susceptibility [14] . It has also been used in our calculation of up 
and down quark mass through the analysis of pion mass and decay constant [15| • Extension 
of these works to the 2+1-flavor case is in progress. 

By using the value of Z^ s we quoted in this article, we are planning to determine the 
up and down quark mass m U( i and the strange quark mass m s from the analysis of the 
meson masses m 2 and m 2 K and the decay constants /„- and fx in the Nj = 2 + 1 dynamical 
simulation. A preliminary results from this project was reported in 161 ]. 
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APPENDIX A: PERTURB ATIVE MATCHING 



In this appendix, we present the details of our matching procedure. 
The matching of an operator O between the MS scheme and the RI/MOM scheme is 
written as 



G MS (/i) 



w 



RI/MOM 
O 



RI/MOM (/io)) 



(A.I) 



where the conversion factor Wq(/i) from a given scheme X to the so-called Renormalization 
Group Invariant (RGI) value is written as 



70 



1 + (7i - A70) 



os(/x) 

47T 



as(jj) 

47T 



+| ((71 - A7o) 2 + 7 2 + - hii - /9 2 7o) ( 
+ (|(7i - A7o) 3 + I (7i - A7o)(72 + Pilo - P1I1 ~ A70) 

+|(73 - Pflo + 2AA7o - klo + Rli ~ Mi - A72)) 



4-7T 



,(A.2) 



to the four-loop order in terms of the strong coupling constant 0:5 (//). (For the coupling 
constant, we always use the MS scheme.) The coefficients A and 7, are given in terms of 
the coefficients of the /5-function /3(as) and the anomalous dimension 70(0:5) 



p ~ " A (4^ " ^ 2 j^r " 



70 



7 <? 47T 7 ° UJ 7 ° UJ 10 UtT 



(A.3) 
(A.4) 



as A = A/A and 7* = 7$ /A- 
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The /3-function is specified by 

2 



A) 



11 - -JV/, 
3 

102 - f AT /5 

2857 _ 5033 325 
2 18 7 54 



149753 



6 



+ 3564C 3 - 



1078361 6508 



162 



27 



+ 



/ 50065 6472 



I" 



162 



+ 



Ca ) Nf + 



1093 Ar3 

N*, 

729 J 



(A.5) 
(A.6) 
(A.7) 

(A.8) 



with ( 3 = 1.2020569. The running coupling constant is then obtained as 



47T 



1 - 



0o In /A 2 

+ # 

/3 4 (ln/x 2 /A 2 ) 2 

A 3 

/3 6 (ln// 2 /A 2 ) 3 



A ln(ln/x 2 /A 2 ) 
/f ln/x 2 /A 2 



(lnln/i 2 /A 2 ) 2 -ln (ln^/A 2 ) + 



P2P0 
Pi 



2\2 



(lnln/i 2 /A 2 ) -^(lnln/i 2 /A 2 ) 



3/?o& 
PI 



lnln/i 2 /A 2 + - - 



2A 3 



(A.9) 



In our work, we chose A = 245 MeV for both 2 and 2+1-flavor analysis. 

The renormalization of the scalar bilinear operator O = S = qq is an inverse of the 
mass renormalization. The anomalous dimension thus has a relation 7 m = —7s- At the 



lowest order, 
four- loop order in the RI/MOM scheme 

A 1 ) 



4 for any scheme. Higher order coefficients are calculated in 



30| to the 



X 



T (2) 
' m 



52 

126 Nt, 

9 

20911 3344 



(3 



18386 128 



27 



+ ^rCs Nf + 



9 



928 



N 



(A.10) 
(A.ll) 



7-, 



(3) 



300665987 15000871 _ 6160 „ 

<3 + -^Cs + 



648 108 
/ 670948 6416 
V 243 27~ 



Cs N - 



3 

18832 
729 



7535473 627127, 4160 \ Ar 
<3 + — Cs N f 



108 



54 



Nj, 



(A.12) 
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and in the MS scheme 

^ = lata -(^ + f(,)y f -^N), (A.14) 

, 3l 4603055 X 35680 . , /91723 34 192 ooM , 18400 . \ Ar 

7< 3) = + — Cs - 8800 C5 - + — Cs - 880C 4 - — Cs J 

/5242 800 160 \ Ar2 /332 64 \ , . , 

where Cs = 1.202057, (4 = vr 4 /90 and Cs = 1.036928. 

For the quark field renormalization (O — q), the lowest order coefficient vanishes in the 
andau gauge for any scheme. The higher order coefficients for the RI/MOM scheme are 



30] 



= (A.16) 



3 3 



m 43477 607 A /3674 \ Ar 80 Ar2 /a 

^ 54714743 7004309 , 15846715 , Z4659455 637413, nort> \ Ar 

r } = Cs + Cs Cs + 830C 5 N f 

648 162 s 1296 \ 324 162 J 

( 166269 „,. \ Ar2 688 Ar3 
+ ( ?T7o " 64 Cs iVj - — iV), (A.18) 



\ 243 

while the MS coefficients are 

1? = y - (A.19) 

(2) 20729 79 . 550 Ar 20 Ar2 

Tn = Cs N f + — N f . (A.20) 

lq 36 2 S3 9 t 27 7 ' V ; 

(3) 2109389 565939 2607 761525 /162103 2291 79 160 \ 

% - ^Cs + — C 4 - y^Cs - + ^Cs + yC4 + — Cs j N f 

/3853 160 \ Ar2 140 Ar3 

+ + Cs Nf + Nf. (A.21) 

V 81 9 / 243 

Coefficients for the tensor current O = T = qo^q are known to three-loop [sl,[3l|. Besides 

the common value 7^? = 4/3, higher order coefficients for the RI/MOM scheme are 

(1) 362 52 Ar 

IF = — ~ 27#/, (A.22) 

(2) 159607 130 72 / 17426 256 \ AT 928 Ar2 /A , 
It = Cs + + Cs N f + Nf, A.23 

81 27 S3 V 81 27 s3 / / 243 f ' V ; 
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and those for the MS scheme are 



(1) 362 52 Ar , kni s 
1? = — - zjNf, (A.24) 

(2) 52555 928 /5240 160 \ 4 2 
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